Burgulence and Alfven waves heating mechanism of solar corona 
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Heating of magnetized turbulent plasma is calculated in the framework of Burgers turbulence 
[A.M. Polyakov, Phys. Rev. E. 52, 6183 (1995)]. Explicit formula for the energy flux of Alfven 
waves along the magnetic field is presented. The Alfven waves are considered as intermediary 
between the turbulent energy and the heat. The derived results are related to a wave channel 
of heating of the solar corona. If we incorporate amplification of Alfven waves by shear flow the 
suggested model of heating can be applied to analysis of the missing viscosity of accretion discs 
and to reveal why the quasars are the most powerful sources of light in the universe. We suppose 
that the Langevin-Burgers approach to turbulence we have applied in the current work can be also 
helpful for other systems where we have intensive interaction between a stochastic turbulent system 
and waves and can be used in many multidisciplinary researches in hydrodynamics and MHD. 

PACS numbers: 



I. INTRODUCTION 

For more than 60 years we are still facing the perplex- 
ing astrophysical problem why the temperature of the 
solar corona is two orders of magnitude higher than the 
temperature of the photosphere; for review and refer- 
ences see for example Ref. [3, sec. 5.1]. The purpose of 
the current work is to investigate the wave mechanism 
of heating of the solar corona, according to which the 
energy is being transported by magneto-hydrodynamical 
waves, generated in the stochastic convective region. In 
the framework of this scenario the famous correlation be- 
tween the solar activity and the coronal emission in the 
X-ray range can be easily explained. The novelty in the 
model is the incorporation of the Burgers approach to 
the turbulence in the problem for the calculation of the 
spectral density of the Alfven waves. In the framework 
of this model we derived an explicit formulas for: 1) the 
spectral density of the MHD waves, 2) the total energy 
density and 3) the dissipated power per unit volume. In 
perspective the results derived here can be generalized 
for a realistic model of the turbulence and the distribu- 
tion of density, temperature and magnetic field. In the 
current work the model task for calculation of the wave 
power, emitted by a turbulent half-space and transmit- 
ted in a non-moving magnetized plasma, is considered 
as an illustrative example. The homogeneous magnetic 
field is oriented perpendicularly to the interface of both 
half-spaces. 



II. BURGULENCE 

The Burgers approach to the turbulence gives an op- 
portunity for approximate treatment of plenty of prob- 
lems; as a tutorial lecture on Burgulence and a source of 
main references we recommend Ref. [5j]. This approach 
to the turbulence has been introduced in the astrophysics 



by Zeldovichi, but didn't become very popular. Yet in 
the last decade the incorporation of the quantum field 
theories gave incitement to the theory of turbulence. 
Kolmogorov^ power laws have been derived and grad- 
ually the Burgers approach converted from a sophisti- 
cated high energy physics problem^ to a standard tool 
for investigation of the turbulence over different physi- 
cal phenomena. The concept of random driver in the 
special case of formation of spicules and the problem of 
heating of space plasma in general was considered in the 
recent works Ref. [l5| : here we would like to emphasize 
on the similarity between propagation of Alfven waves 
in an infinite magnetized plasma and kink waves in solar 
spicules^; for other papers devoted on the subject see 
Ref. [3] • There exist also different ways of generation of 
waves by turbulence where the influence of the random 
driver is negligible 12 . 

In the next section we will consider the application of 
the Burgers approach to the theory of creation of Alfven 
waves from the turbulence and the stochastic granula- 
tion. In the (3 ~ 1 region of the solar atmosphere Alfven 
waves can also be generated by resonant conversion of 
acoustic oscillations^!. 



III. WAVE AMPLITUDES AND SPECTRAL 
DENSITIES 

According to the Burgers idea£ the influence of the 
turbulent vortices on the fluid can be modeled by intro- 
duction of a random volume density of an external force 
F(t, r) on the right hand side of the Navier-Stokes equa- 
tion 

p(<9 t +V- V)V = -Vp + 7?AV+j xB+F. (1) 

Analyzing low frequency and long wavelength phenom- 
ena one can approximately suppose the random force to 
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be a white noise with a <5-function correlators^ 

(F(t 1 ,r 1 )F(i 2 ,r 2 )}=fp 2 5(i 1 -i 2 )5(r 1 -r 2 )l, (2) 

where T is the Burgers parameter. For the other notions 
we use standard notations: p is the pressure, r\ is the 
viscosity, B is the magnetic field. 
The current j is given by Ohm's law 



j = o-E', E' = E + VxB, 



(3) 



where a is the electrical conductivity and E' expresses 
the effective electric field acting on the fluid. 

In quasi-stationary approximation, when the <9fE term 
is negligible, the evolution of the magnetic field in a 
plasma with constant resistivity is 

d t B = rot (V x B) - ^ m rotrotB. (4) 

We consider an incompressible fluid 

divV = (5) 

with wave amplitudes for the velocity 

V = (V x ,V y ,V z ), V^V A (6) 

and the magnetic field 

B = B + B', B'(t,v) = (B' x ,B' y ,B' z ), B' <C B (7) 

being small compared to the Alfven speed and the ex- 
ternal magnetic field. The z-axes is chosen in a vertical 
direction along the constant magnetic field 



B = (0,0, S ) = B e z 



(8) 



Analogously for the pressure we also suppose small devi- 
ations from a constant value po 



P = Pa+P 



(9) 



The linearized system of magnetohydrodynamic (MHD) 
equations then reads 

V»' F Bn ( dzB ' x ~ dxB z 
d t V=-^+ - + ^°- d z B' y - d y B' z |+ ^AV, 
P P PoP\ o 



a t B'=s a 2 v + ^ m AB', 

divB' = 0, divV=0, 



where 



_ r\ 

Vk = ~ 

P 



(10) 



(11) 



are respectively the magnetic and kinematic viscosities. 
Here we would like to emphasize that the incompress- 
ibility condition divV = is applicable in the convec- 
tive zone of the sun, where the sound speed significantly 
exceeds the Alfven speed c s 3> Va- In this region the 
convective circulation is treated as the influence of the 



turbulence in the Burgers approach. The excited Alfven 
waves then propagate along the magnetic filed lines and 
reach the corona where the direction of the inequality 
could be opposite. 

Let us consider the evolution of the amplitudes of 
standing plane waves 

V(t,r) = V A v k (t)sm(k-r), (12) 

B'(t,r) =B bk(*) cos(k-r), (13) 
F(t,r) =pV A f k (t) sin(k-r), 

p = pVhP k cos(k • r), (14) 

(£(ti)fc(t2)) = T6(h - t 2 )S pM l. (15) 

This is a technical tool for evaluation of the average 
energy density used in many textbooks on statistical 
physics. Due to the negligible dependance on the bound- 
ary conditions usage of running waves would lead to the 
same result. 

Substitution of this ansatz in Eq. (|10j) gives a separa- 
tion of the variables and for the wave amplitudes we have 
a system of ordinary differential equations 



d,c' = //k + ./' - I a j k-hy - k y b z ] - i/fefcV. (Hi) 

d t h = V A k z v - v m k 2 b, (17) 
k-b = 0, (18) 
k-w = 0, (19) 

where for the sake of brevity the wave- vector indices are 
omitted. Time differentiation of the incompressibility 
equation Eq. (fT9|) and substitution in Eq. (fTB)) gives the 
explicit form of the pressure 



kf 

P = — p Vao z 



(20) 



Further the substitution of the pressure from Eq. (f2"0")) in 
Eq. (fl6|) gives 



where 



v = —VAk z h — v\Jt 2 v + i^, 
b = V A k z v - v m k 2 h, 



fx = III 



(21) 



(22) 



is the transverse projection of the external force with 
respect to the wave- vector k by the polarization operator 
Hx 



n 



k® k 



k 2 y +kl 
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k l+ k l 
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k l+ k l 






fe2 


k' 1 


k' 2 



(23) 



This method for elimination of the pressure is used in 
Ref. [ll[ dedicated on amplification of Alfven waves in 
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the presence of shear flow. A continuation of this research 
in the axial symmetric case is presented in Ref. (l3| |. 

Let us investigate the eigen-modes oc exp(— itut) dis- 
carding for a moment the influence of the external force 
f . For the case of small dissipation the secular equations 
for all the components are the same 



-iu> + v\Jt 2 —Vxk z 



= 0, 



with eigen-values 



(24) 



(25) 



where ua = Vpjt z is the frequency of the Alfven waves 
and 



7 = vk 2 , v = i> n 



(26) 



is the attenuation coefficient. 

In regime of negligible magnetic viscosity the MHD 
system Eq. (1211) turns into the second order ordinary dif- 
ferential equation 

v = v^k 2 v — lo\v + fi. (27) 

If we introduce the ^-component of the displacement 



x(t) 



v x {t')dt' 



(28) 



we obtain an effective oscillator equation under an exter- 
nal force 

x = v-^k 2 x — lo\x + fj_, (29) 

which in case of negligible kinematic viscosity has an ef- 
fective energy 



£ = \ if + <* 2 ) 



(30) 



This is the non-perturbed by the friction and the weak 
random noise energy density. A harmonic oscillator un- 
der an external force is a well-known mechanical problem, 
see Ref. fl, Chap. V, Eq. (22,11). When the external 
force is a white noise a simple averaging gives that the 
energy is a linear function of time; in other words the 
random noise gives a constant power to the oscillator 

d t (S) = d t + \ufo*} = d t (v 2 ) = E, (31) 



where 



r 



(32) 



is the Burgers parameter in the correlator for the plane- 
wave amplitudes of the transverse projection of the ran- 
dom force in Eq. (fT5|) 



(33) 



C,±(*l)fk,±(* 2 )) - r<*(ti - h)S pM t; 



the technical details of this simple derivation will be omit- 
ted. Here we suppose that the wave energy does not 
increase significantly over one period 



(34) 



and the brackets stand for a period and noise averaging. 
In static regime the condition for dynamic equilibrium 
requires that the power received from the noise Eq. (|31|) 
has to be equal to the dissipated power 



d t (£) = - 1 (£) 



(35) 



and this gives the dimensionless static averaged energy 
of the oscillator 



r 



(36) 



Now we are going to use this model example for deriva- 
tion of the spectral density of Alfven waves. For the vol- 
ume density of the average wave energy we have 



(E) = 



(B' 2 ) (V 2 ) 



2fi 2p 
This result may be rewritten as 



(37) 



where we have preformed a volume averaging 

(sm 2 (k.r)) = i, 
and used the definitions for the magnetic pressure 



and Alfven speed 



^3. 

2fl 



YL = *L. 

2p 2p ' 



(39) 



(40) 



(41) 



In Eq. (|38p we can substitute the static squared am- 
plitude (v 2 ) s t from Eq. ([36| and the attenuation coef- 
ficient from Eq. (f2"6"|) . Taking into account the y and 
z-components of the Alfven wave velocities results in an 
additional triplication of the energy. Then for the static 
value of the density of the wave energy we finally derive 



E k = (E) st 



3Tp E 



4 vk 2 



(42) 



This spectral density is the main detail for the statistical 
analysis made in the next section. 
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IV. ENERGY DENSITY, HEATING RATE AND 
ENERGY FLUX 

In order to calculate the total static energy density 
we have to perform a summation of the spectral density 
Eq. ([42]) over all wave- vectors 



d 3 k 

(2^ 



(43) 



The wave- vectors' cut-off is determined by the condition 
for the existence of Alfven waves, i.e. by the region where 
the Alfven frequency wa equals the attenuation coeffi- 
cient 7 



V A k z = uk 2 



(44) 



This equation sets the natural cut-off k c for the vertical 
component of the wave- vector k z , parallel to the constant 
magnetic field B 



kr . 



(45) 



Taking into consideration that the Alfven waves are 
spreading axially symmetric as they follow the magnetic 
field lines we may rewrite the wave- vector as 



k 2 = kl 



(46) 



which fixes the maximal value of the plane components 
of the wave-vect< 
ponent cut-off k c 



of the wave-vector k p as a function of the vertical com 



£,(max) 



\fkjk c 



(47) 



These cut-offs will be used as boundary conditions in the 
calculations of the total wave energy, integrated over all 
possible wave- vectors in Eq. (|43p . and the corresponding 
wave power. In this way considering the axial symmetry 
we obtain 



{2n)\ A.=i Jk 



1 



d(nk 2 )dk Zl (48) 



where we have taken into account that due to infrared 
divergencies the vertical component of the wave-vector 
should also have a minimal value, determined by the typ- 
ical size R of the given magnetic field loop. Then, having 
in mind Eq. (|32p , the final expression for the total energy 
becomes 



STRp B 
8tt 2 ^ 3 



In 



V A R 



(49) 



Hereby, for the Alfven waves energy flux along the mag- 
netic field we find 



1 _ 3TV A Rp E 

2 totVA ~ 16^ 



-(¥)-' 



(50) 



Multiplication of the spectral density with decay 7 gives 
the heating rate and summation in the wave-vectors' 
space gives the dissipated by the Alfven waves per unit 
volume power 



Q = V / E klk 



d 3 k 



TV A 



(51) 



Perhaps magnetic field dependance of the heating rate 
can explain the correlation between the magnetic flux 
and X-ray luminosity^. In order to apply this result we 
need to know the correlation between the magnetic field 
and the radius of the active region. 

The parameter with dimension length R which we in- 
troduced due to the infrared divergencies participate in 
the final formulas for the total energy density Eq. (j4"9"|) 
and energy flux Eq. ([50]) . It seems very plausible if this 
detail could be helpful to reveal the scaling correlations 
of the length of coronal loops and other observable quan- 
tities. It is premature, however, to start a detailed dis- 
cussion based on a model example. The purpose of the 
present work is to demonstrate that the Burgers approach 
should be used for realistic computer simulations in the 
future. 



V. DISCUSSION AND CONCLUSIONS 

Up to now, the coronal heating problem has not a sat- 
isfactory explanation not because of a lack of interest, 
but rather due to difficulties concerning the simultaneous 
obtaining of reliable observational data for all processes 
running on the solar surface. For determination of the 
parameters of each model as a rule we are forced to use 
indirect means. The model, which we now present on the 
arena suffers the same disadvantage. For its confirmation 
it is necessary to know: 1) the distribution of magnetic 
field near the sunspots, 2) velocity-velocity correlator in 
horizontal direction or, which is equivalent, the correla- 
tor of the small stochastic perpendicular components of 
the magnetic field, etc. In short, apart from a realistic 
model for the turbulence we also need a realistic model 
for all other variables. This naturally supposes a broad 
interdisciplinary collaboration, which involves the whole 
spectrum of the solar physics, from the statistical MHD 
to the X-ray luminosity of the solar corona. The confir- 
mation of the suggested model requires realistic Monte 
Carlo simulations, for which the analytical calculations 
done in the current work Eq. 1(50]). Eq. !j4"9"]). Eq. f5l]) 
are only test examples indispensable for a more realistic 
treatment of the problem. We propose an oversimplified, 
and therefore solvable, model, which realistic generaliza- 
tion will hopefully lead to an adequate theory for the 
solar corona heating. 

A natural continuation of the current investigation is 
to take into account the influence of the shear flow over 
the magnetized turbulent plasma. Preliminary numeri- 
cal calculations on the behavior of MHD waves in a shear 
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flow show a significant amplification of the slow magneto- 
sonic mode. The waves energy amplification is actually a 
transformation of the shear flow energy into wave energy 
which at the end dissipates into heat. This dissipation 
originates an effective viscosity in the shear flow and the 
current scenario for heating of the solar corona may also 
happen to be the theory for the missing viscosity in ac- 
cretion disks. Understanding the origin of an additional 
viscosity in a magnetized turbulent plasma is rather im- 
portant for the proper explanation of the shining mech- 
anism of quasars, as well as for the angular momentum 
transport during formation of compact astrophysical ob- 
jects. 

The theory for the coronal heating takes into consid- 
eration the influence of the stochastic random noise of 
the turbulence and the granulation over the evolution of 
Alfven waves. A similar mathematical scheme is used in 
the condensed matter physics for calculation of the fluc- 



tuational superconductivity 4 . In the physics of supercon- 
ductors the wave field is the effective Ginzburg-Landau 
wave function and the noise is produced by Langevin 
thermal fluctuations, which come from the inhomoge- 
neous term in the TDGL equation. Naturally, the solid 
state scheme is meticulously developed and the final ana- 
lytical results can be used for fitting of repeatable exper- 
imental data from current-voltage characteristics. The 
present work represents a realization of a similar idea in 
the new astrophysical plasma areal. 
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